CERN-TH/98-55 
IFT/3/98 



Higher-order QED corrections to e + e — > 1/1/7 at LEP2 



A. Jacholkowska 1 , J. Kalinowski 2 and Z. Was 3,4 

1 Laboratoire de VAccelerateur Lmeaire,CNRS-IN2P3, 914050 ORSAY, France 
2 Institute of Theoretical Physics, Hoza 69, 00681 Warsaw, Poland 
3 CERN, Theory Division, CH-1211 Geneva 23, Switzerland 
* Institute of Nuclear Physics, ul. Kawiory 26a, Cracow, Poland 



Abstract 

The process e + e~ — > vv + 7 with a distinctive "photon-plus-missing-energy" 
signal can serve as an important tool to search for new physics at LEP2. It can be 
exploited to measure the WWj coupling, or to search for weakly interacting and 
invisible (s)particles. For meaningful comparisons of experimental data with theo- 
retical predictions, higher-order QED corrections due to multiphoton emission must 
be taken into account. In the present paper we explain how the WWj coupling has 
been incorporated into the KORALZ Monte Carlo program, which can now be used 
for simulations of e + e~ — > vv + 7 events with higher-order QED corrections. The 
strategy of how to perform an experimental analysis in the presence of experimental 
cuts is proposed. The question of systematic uncertainties is addressed and some 
numerical results for the phenomenologically interesting case of anomalous WW'j 
couplings are also given. 
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1. Introduction 



One of the main goals of experiments at high energy consists in comparing new and 
so far unexplored measurable quantities with the Standard Model (SM) predictions and 
thus providing its important tests in a new experimental domain. On the other hand, 
any discrepancy can be a sign of new physics. From the practical point of view it is thus 
important to confront new experimental data with theoretical predictions for observables 
which are: (i) sensitive to expected new physics, (ii) under control of systematic un- 
certainties for the data and theoretical predictions, and (iii) have high statistics. To 
optimize the experimental analysis and to suppress SM backgrounds, complicated pat- 
terns of kinematical cuts are often introduced. One should also bear in mind that the 
detector sensitivity varies over detector acceptance in a complicated manner. That is 
why high-precision physics Monte Carlo programs, including complete Standard Model 
predictions as well as contributions from possible 'new physics', are useful and in some 
cases even mandatory for such projects 0. 

The measurement of the nature of three-gauge-boson couplings has been identified as 
one of the most important topics in the LEP2 scientific program 0. The most discussed 
e + e~ — > WW process in this context has a theoretical drawback as it involves two WW Z 
and WW 7 couplings intertwined, leading inevitably to a model-dependent analysis. In 
our study we will concentrate on the phenomenology of a single WWj coupling. In |2], |3| 
it was argued that the 

e + e _ — > vv^ (1) 

channel is a good candidate for such measurements, because the process is sensitive only to 
this coupling, is of clean hadron-free signature and has a rather large cross section. Such 
final states have been used at PETRA, SLAC and LEP1 M as a means of measuring 
the number of light neutrinos since, at these low energies, the contribution from WW^f 
coupling is negligible. However, at higher energies the t-channel VF-exchange process 
becomes important, allowing for a direct study of WWj coupling, independently of the 
WW Z coupling. Therefore it opens an exciting possibility since new physics can, and in 
general will, manifest itself in corrections to the gauge-boson sector, in particular in the 
WWj vertex. 

The events with photon(s) plus missing energy in e + e~ collisions might originate also 
from other mechanisms, signalling new physics beyond the Standard Model. For example, 
such final states can be produced in both gravity- and gauge-mediated supersymmetric 
models. The missing energy in these events is caused by weakly interacting supersymmet- 
ric particles, such as gravitinos, neutralinos and/or sneutrinos M. In all such cases the 
Standard Model e + e _ — > z/z/7 events are irreducible background and reliable predictions 
for them are therefore necessary. 

Recently the process e + e~ — > vv^j has been exploited by the ALEPH collaboration 
to derive preliminary limits on anomalous WWj couplings from data collected at 161, 
172 and 183 GeV 0. Although the results for the couplings Ak 7 = 0.05 (stat) and 
A 7 = —0.05 ii*! (stat) are not competitive with the ones expected to be derived (in a 
model-dependent way) from WW final states at these energies, they demonstrated the 
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physics potential of processes with isolated photons. Process (|T|) has also been used in 
searches for supersymmetric signals 0. 

To draw any conclusions on new physics from the experimental data, however, a proper 
treatment of higher-order QED corrections to the process in eq. (0) is required. In the 
present paper we will discuss a Monte Carlo program that can be useful in reaching such a 
goal. As a specific example we consider the vv^j final states as a means to investigate the 
WWj coupling. Our paper is organized as follows. First we will present basic distributions 
obtained from single bremsstrahlung analysis prior to this work and recall that they were 
shown to be sensitive to anomalous WWj couplings. We will also show that they are 
prone to higher-order QED corrections. Then we will present the KORALZ Monte Carlo 
program in which a realistic algorithm of how to match the single-photon matrix element 
for e + e~ — > vv^i (including anomalous WW'y coupling) with higher-order QED corrections 
has been implemented. We will present consistency tests of such an approach. Finally 
we will turn to the question of the choice of physical quantities that are most sensitive to 
anomalous couplings when higher-order QED corrections are switched on. Discussion of 
numerical results will conclude the paper. 

2. WW 7 vertex and higher-order QED corrections 

Expressed in purely phenomenological terms, the effective Lagrangian for the WW'y 
interaction can be written in terms of seven anomalous form factors ||. Assuming C and 
P conservation^ reduces their number to three, which are conventionally denoted by AgJ, 
Ak 7 and A 7 . The effective Lagrangian thus restricted can be written as follows 

(1 + AgJ)(W^W a - W^W aP )AP + (1 + An^WlWpA^ 

(2) 

where A a p = d a Ap — dpA a and W a p = d a W/3 — dpW a . All higher derivative terms can 
be absorbed by the above couplings, provided they are treated as form factors and not 
constants. 

For the photons on-shell, as is the case in the process of interest e + e~ — > vv , y, electro- 
magnetic gauge invariance requires that Agl(q 2 = 0) = 0. Therefore we are left with only 
two unknown form factors. In the Standard Model both couplings Ak 7 and A 7 vanish. 

In the process e + e~ —> ^^7, there are only two kinematic measurable variables at our 
disposal, for example the energy and the polar angle of the observed photon. In order 
to achieve the best sensitivity to anomalous couplings, it was shown in Refs. [|, ||] that 
double-differential distributions of photons measured in the barrel detector are best suited. 
However, for the discussion of our procedure it will be sufficient to refer to one-dimensional 
differential distributions. 

Subject to possible cancellations between various contributions, data on the neutron electric dipolc 
moment constrain these form factors very severely. 
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We use the following set of selection criteria, which have been used by ALEPH: (i) pho- 
ton transverse momentum pt > 0.05-y/i, (ii) photon energy E 1 > 0.1 >/s, (hi) photon 
angular acceptance region | cos(9 7 | < 0.95, and for definiteness we take the centre-of-mass 
energy y/s = 172 GeV. In Fig. [I] the energy, the angular and the transverse momentum 
distributions for a single photon are shown for three cases: genuine O(a) SM results (solid 
lines), SM including multiphoton radiation but with WWj interaction exluded (stars), 
and genuine O(a) including anomalous WW'y coupling (open circles). The maximum 
around E r = 60 GeV in the energy distribution corresponds to the so-called radiative 
return to the Z, i.e. the photon energy for which the invariant mass of the recoiling 
system vv pair is close to the Z mass. 

There are several comments to be made. Comparing solid lines with stars, it is seen 
that photons with energies higher than E r are best suited for the WW'y coupling mea- 
surement, because of a sizeable sensitivity to anomalous contributions. For the energies 
matching radiative return, the significance is very small because of the high Standard 
Model background. For energies below E r the photons are also sensitive to anomalous 
couplings, although the sensitivity decreases systematically with decreasing photon en- 
ergy. Below 20 GeV the sensitivity is negligible even to rather strong couplings, which we 
have taken in the calculations, namely Ak 7 = —10, A 7 = 0. We will use this set of anoma- 
lous couplings throughout the paper as an example. The above conclusions are derived at 
the single-photon emission level. This might be sufficient to estimate the discovery power 
of the appropriate observables. However, it must be completed with a careful analysis 
of higher-order QED corrections because these may change not only the size but also 
the shape of the distributions that are analysed. To illustrate this point, in Fig. [I] we 
include for comparison the results of a multiphoton higher-order Monte Carlo simulation 



performed with the help of the KORALZ program [10], in which the WW'y coupling is 
not included (open circles). For photons of energies higher than that of radiative return, 
the effect of higher-order QED corrections is small; they slightly decrease the SM cross 
section. In that kinematical region the inclusion of higher-order correction does not seem 
to be essential for our choice of the anomalous couplings. The situation is different for 
photon energies below that for radiative return. Contributions from both anomalous cou- 
plings and higher-order QED corrections to the distributions are of similar order and have 
the same sign (for the values of Ak 7 and A 7 chosen). Higher-order corrections vary in 
that region from —30% to about 90%. This result shows that any realistic determination 
of the anomalous WW'y coupling must compare data with theoretical predictions that 
include anomalous terms as well as higher-order QED corrections at the same time. This 
observation is the main motivation for our present work. 

3. Monte Carlo generator 

The most up-to-date published description of the KORALZ 4.02 program for produc- 
tion of pairs of leptons, including radiative QED corrections, is described in detail in 
Ref. [TIJ. Real photon radiation can be generated at different levels of sophistication. 
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Figure 1: Inclusive photon distributions for: (i) genuine order-a distributions in the 
Standard Model - solid line; (ii) order-a including anomalous interactions - open circles; 
and (iii) the SM with multiphoton radiation - no WWj interaction - stars. For the set 
of experimental cuts and anomalous couplings see the text. 



First, the program may be run at the Born level for the e + e _ — > 11 process^. Second, it 
may be run at order a when single bremsstrahlung configurations are generated. Finally, 
the program can be run with pragmatic order-a 2 QED corrections including exclusive 
exponentiation. In this case the number of photons actually generated is not restricted. 

Let us briefly recall the main features of KORALZ existing prior to the present work, 
which are relevant for e + e _ — > vv^i production. The basic design of KORALZ is for the 
s-channel processes, such as muon or tau lepton pair production. The cases of muon 
and tau neutrinos accompanied by a single photon, see diagrams (1) and (2) in Fig. 0, 
have been implemented easily, since the modification is straightforward, consisting of 
changing the appropriate coupling constants only. In the case of electron neutrinos, there 
are additional contributions from i-channel W exchange, diagrams (3) and (4), and from 
the WW j vertex, diagram (5) in Fig. ^. Since in physical gauge the amplitude for the 
diagram (5) with WW-j vertex is finite in soft and collinear limits and negligible at LEP1 



11 1, it was not implemented at the time. Note that by neglecting the WWj vertex 



the generated distributions are formally not gauge- invariant. Diagrams (3) and (4) have 
2 This option is not interesting for I = v since no photon radiation is present. 
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Figure 2: T/ie Feynman diagrams for e + e —>■ vvy. 



been implemented in an approximate way[] as follows. First a series of (weighted) events is 
generated inside KORALZ according to the Born-level s-channel Z-exchange diagrams for 
neutrino pair production (leading to totally unobservable final states), and with higher- 
order QED corrections including additional virtual and/or real photons. In this way final 
states consisting of 0, 1, 2, 3, photons accompanying the neutrino pair are generated 
according to the prescription of exclusive exponentiation of YFS [12 . Full phase space 
is covered by the generation. At this point the complete event with four momenta of 
final-state neutrinos and photons is constructed and stored in the common block. This 
procedure is correct for v^v^n^ final states. For the electron neutrinos, the matrix element 
for the hard process e + e~ — > v^v^ (with Z exchange only) has to be replaced by the matrix 
element for e + e~ — *> v e v e (with Z and W exchanges). Consider first 0(a) contributions 
coming from Z (1+2) and W (3+4) diagrams. Thanks to the factorization properties of 
the QED, the matrix element squared |A4o| 2 for the sum of four diagrams (1+2+3+4) is 
approximated as a product of the matrix element squared I.M12I 2 for two diagrams (1+2) 
and the weight factor F w : 



IM 



1 w 



I approx J 2 

\M Z + M 



\M 



12 



Wi 



\v\ 



(3) 



where M. z and Aiw denote amplitudes for v e v e pair production via Z- and W-exchange 
diagrams, respectively, with no photons in the final state. In the calculation of the factor 
Fw four momenta of v e and v e are used in a straightforward way, but for the incoming 
electrons and positrons effective beam momenta are calculated from momenta of the real 
beams and real photon emitted, according to the prescription given in Refs. 10, 13 



3 It is thus important to realize that the program quality for v e v e l final state was (and remains) 
inferior to that for other final states. 
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Note that the factor Fw given in eq. (H) is gauge- invariant. The standard final rejection 
of KORALZ creates the sample of unweighted events, which is passed to the program 
user. The factor F\y is included in that final rejection weight. 

Let us stress that in this way the effect of diagrams (3) and (4) with t-channel W 
exchange, as well as their interference with the diagrams (1) and (2), has been introduced 
not only for case of a single-photon radiation, but also for multiphoton configurations. 
This approximation, as well as the whole program, has been heavily tested at LEP1 
energies both for the single bremsstrahlung mode and the multiple photon case. Extensive 
comparisons with the results from another program [fH|| , including a complete matrix 



element for single and double bremsstrahlung, have been performed, finding an excellent 
agreement and thus proving both the validity of the approximation and the correctness 



of the code, see Refs. [TO, |13|, 15] for more details. During the last seven years, i. e. since 



the vv pair production process has been implemented into KORALZ for the first time 



151 , the program was modified several times, and in general its applicability for s-channel 



processes was assured also at LEP2 energies, see [jTS . 

For the case v e v e ^^ however, the WW'j interaction (and the question of gauge invari- 
ance) can no longer be ignored, and the approximate treatment of t-channel W exchange 
requires refinements to be precise enough at LEP2 energies. 

We now turn to the discussion of how these problems are solved in the present ver- 
sion 4.04, which is the modification of the version 4.02 of the KORALZ program. The 
modification consists of introducing the new weight defined as a ratio of matrix element 
squared for e + e~ — > vv^j with and without the WW^ vertex:^ 



= \M Q + M WWl \ 2 
l-Mol 2 



The term Aiww-y contains contributions from diagram (5) with the WW^f vertex, in- 
cluding anomalous couplings. The matrix element is taken from the numerical program 
developed in Ref. H]. Note that the term in the denominator |.Mo| 2 is gauge dependent. 
If in the present version KORALZ 4.04 the events were generated according to the matrix 
element |.Mo| 2 , the gauge dependence would cancel out in the product of eqs. (^) and 
However, this is not the case, since they are generated according to the approximate 
\M.'q >ptox \ 2 as explained above. Nevertheless, as we will discuss in the next chapter, see 
Fig. the numerical error is small and the final distributions are affected by this ap- 
proximation in a small and controlled manner only. The implementation of the above 
weight Wt proceeds as follows. For every event, the complete kinematic configuration is 
taken. In cases when there is just one photon, the formula ([D can be used in a direct way. 
However, for events with multiple photons generated, there are many possibilities for a 
choice of a photon radiated off the W boson. Since the QED bremsstrahlung of high-p^ 

4 In some regions of phase space, due to destructive interference of contributions from the s- and t- 
channel boson exchanges, the Standard Model predictions can be very small. There the weight defined in 
eq. (|J) can in principle tend to infinity, making the simulation numerically non-convergent. This technical 
obstacle, not relevant to our discussion, is overcome by replacing the SM matrix element by a different 
one and obtaining the SM predictions also in the form of a weighted event sample. 
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Figure 3: The ratio of distributions generated by the KORALZ 4.04 and the program 
used in Ref. ||. Both programs are run with the same input parameters and kinematical 
cuts. Anomalous WWj couplings are included. 



photons from electrons is strongly suppressed, and the WW'y contribution is important 
for energetic photons, in our algorithm we calculate^ the weight factor Wt in eq. (H) only 
for the photon of highest pt- All remaining photons are then incorporated into effective 
beams according to the following recipe: 

a) Except for the photon of highest pr, all remaining photon's four-momenta are sub- 
tracted from the four-momentum of the beam of the same hemisphere. In this way we 
obtain effective electron beams, which are not necessarily on-mass-shell. 

b) We boost all four-momenta of effective beams, neutrinos as well as the highest-jo^ pho- 
ton to the rest frame of effective beams. 

c) We reset the four-momenta of the effective electron beams to be on-mass-shell in such 
a way that their back-to-back direction and energy-momentum conservation is preserved. 

At this point we have constructed reduced on-mass-shell kinematics for e + e~ — > vv^j 
configuration for which the weight Wp can be calculated according to eq. (^). For practical 
reasons explained at the end of the next chapter, we do not include the Wt weight into 

5 At this point we assume that the WWj couplings do not affect the matrix elements in soft or collinear 
photon limits, or virtual corrections. In short, it means that the whole structure of QED remains intact 
and that the sole effect of new couplings is a modification of the amplitude for configurations with photons 
of high energy and px- 
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KORALZ internal rejection loop. Instead, generated sample follows the same distributions 
as in previous versions of the program. The WW'j contribution manifests itself only 
through the weight Wt which is provided for every event. 

4. Numerical tests 

The reduction procedure presented above is not free of systematic uncertainties for 
WW^j contributions to multiphoton production processes, except at the single photon 
level. These uncertainties, though, are beyond the leading-log QED level as well as beyond 
soft-photon corrections. Only genuine next-to-leading-log corrections to contributions 
from WW'j interactions are out of control. On the other hand, one benefits from full 
coverage of the bremsstrahlung photons phase space provided by the host Monte Carlo 
as well as its predictions for Standard Model without any deterioration in quality. 

As a first step in testing the implementation of the WW'j vertex (with or without 
anomalous couplings) and diagrams including t-channel jy-exchange into KORALZ, we 
performed studies at a single-photon level. In this case results from KORALZ should be 
the same as from direct Monte Carlo integration (using weighted event sample), as e.g. 
in the original approach of ||. Any discrepancy would indicate either technical problems 
in any of the two programs, or breakdown of our approach. The results are presented in 
Fig. |3|, where the ratios of differential cross sections obtained from two programs are shown. 
The results are consistent over the entire phase space^. One can see that distributions 
differ by at most a few percent. This also shows that the numerical error due to breaking 
the gauge invariance by our procedure is small. 

In the second step we verify that the reduction procedure to effective beams as ex- 
plained in the previous section indeed works. For that purpose we compare at a given cms 
(center-of-mass system) energy the results from KORALZ with anomalous contributions 
and higher-order corrections included as described above, with another approach (called 
EMU), where convolution of the WW'j amplitude with higher-order corrections is per- 
formed as explained below. In the EMU approach, we first generate the spectrum of M e g 
from the multiphoton run of KORALZ at cms energy without any WW'j interaction. The 
M e ff is the invariant mass of the vv^p-max system. For practical reasons we divide the M e jj 
spectrum into discrete bins. In the second step, for each bin i containing iVj events we 
generate an equal number of events, this time, however, at "an effective cms energy" equal 
to the centre of the bin M e jj with the help of a single bremsstrahlung mode of KORALZ 
operation and contributions from the WW'j interaction included (which we have shown 
before to work properly). In this case additional photons are not present, but we emulate 
their presence by a boost of all final-state particles from the "effective cms energy" M e g 
to the overall cms frame as if there was (just one) additional photon collinear to one of 
the beams. 

As can be seen in Fig. [|, the two plots differ in the physically interesting region|] by 

6 For large E 1 or p7p the cross section is very small. Large statistical fluctuations in the ratios of the 
Monte Carlo results are thus expected. 

7 Physically interesting region consists in this case of photons of pr > 20 GeV and of energies different 
from radiative return to the Z . 
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Figure 4: Comparison of differential distributions obtained from the new version of KO- 
RALZ 4.04 including anomalous interactions (Standard), and the EMU algorithm (trian- 
gles), as explained in the text. 



only up to a few per cent. This indicates that the proposed algorithm is constructed 
properly and also that the bulk of the higher-order corrections is of leading-log type. 

One should bear in mind that such tests are indicative only. To answer the question of 
how good such an algorithm is, would require careful comparisons of the generated distri- 
butions with those from exact matrix-element calculations including anomalous couplings 
of order at least a 2 . This is difficult and not necessary at the moment, as a precision on 
higher-order corrections to anomalous contributions of order 10% is sufficient. 

In practical applications it is convenient to use the method of correlated samples. 
The WWj vertex is linear in Ak 7 and A 7 , the two independent anomalous couplings. 
Therefore the matrix element squared, and consequently all differential distributions, are 
bilinear in the variables Ak 7 and A 7 . It is thus convenient to calculate, for every generated 
point in phase-space, the value of the matrix element at six fixed and conveniently chosen 
combinations of Ak 7 and A 7 values^. From this, the matrix element for any other value 
of anomalous couplings can be easily recalculated. This property extends to differential 
distributions as well. It also remains valid when higher-order leading-log QED corrections 
as explained above are included. 

8 In KORALZ numerical values of (Ak 7 , A 7 ) actually implemented are: (0, 0), (-10, 0), (10, 0), (0, 
-10), (0, 10), (10, 10). 
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Note that KORALZ generates unweighted events for the Standard Model with ad- 
ditional weights corresponding to different modifications of the matrix element due to 
WWj couplings. This not only saves computing time, as it is sufficient to perform detec- 
tor simulation only once, but also the standard and non-standard model samples obtained 



this way are statistically correlated, as was explained in 17[| . This significantly improves 
the convergence of the predictions since the statistical error affects only the difference of 
the two distributions. 

After these modifications, KORALZ version 4.04 (or higher) can be used0 for our 
purpose with any combination of flags for the initial state. Every event generated as 
weight-one event, corresponding to the Standard Model with WWj vertex excluded^ 
is accompanied now by the weight record, allowing for a calculation of weights for any 
numerical value of WW'y coupling, standard or anomalous. This can be especially useful 
in performing final fits to the data in the presence of experimental cuts. It can also be 
used for Standard Model background calculations for other processes with a single photon 
plus missing energy in the final state. 

Finally, let us note that the question of systematic errors of the genuine s-channel 
process for non-electron neutrino pair production was not addressed. It calls for a separate 
discussion including tests/comparisons with, for instance, the matrix element for the 
production of a neutrino pair with up to three real protons. 

5. Results and Conclusions 

Let us now turn to the presentation of the numerical results for anomalous couplings. 
In Fig. |5] we compare the photon spectra normalized to the SM lowest-order calcula- 
tions for three cases: (i) anomalous couplings are included at the lowest order (circles), 

(ii) higher-order QED corrections but no anomalous couplings are included (stars), and 

(iii) both anomalous couplings and higher-order corrections are included (triangles). The 
difference of physical content between this plot and the one of Fig. [I] consists primar- 
ily in presenting the results for the case where higher-order radiative corrections and 
anomalous contributions are included simultaneously. Although the QED higher-order 
bremsstrahlung corrections are sizeable, and varying from —30% to 100%, they do not 
wash out the effect of the anomalous contribution. What is more interesting, in the region 
of high-energy photons the presence of higher-order effects tends to slightly enhance the 
sensitivity to anomalous couplings. On the other hand, for photon energies below that 
of radiative return, the sensitivity deteriorates significantly. We also observe that the 
combined effect of anomalous couplings and radiative corrections is not just a simple sum 
of the two effects. 

We have thus shown that the interpretation of experimental data for anomalous WW'y 
couplings in the e + e~ — > vv^j channel requires theoretical predictions in the form of 
a Monte Carlo simulation program, which includes anomalous contributions as well as 

9 The version 4.04 can be obtained upon request from Zbigniew.Was@cern.ch. It includes also options 
described in Q. 

10 Exactly as in the versions of the program prior to the modifications presented here. 
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Figure 5: The ratios of E 1 and cos 6* 7 distributions obtained from: (i) order-a distributions 
including anomalous couplings - open circles, (ii) multiple photon radiation, no anomalous 
couplings - stars, (hi) multiple photon radiation and anomalous couplings included - 
triangles. In all cases distributions are divided by the appropriate ones obtained from 
the single-photon generation in the Standard Model (no anomalous couplings). All input 
parameters as explained in the text. 



multiple bremsstrahlung higher-order QED corrections. The current version 4.04 of the 
KORALZ program, with modifications presented in this paper, can now be used to achieve 
this goal. It can also be used for Standard Model background calculations for other 
interesting processes, such as gravitino or neutralino/sneutrino production processes at 
LEP2. 
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